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Abstract 

The relic gravitational waves are the cleanest probe of the violent times in the very early history of the Universe. They are 
expected to leave signatures in the observed cosmic microwave background anisotropies. We significantly improved our previous 
analysis [l[ of the 5-year WMAP TT and TE data at lower multipoles £. This more general analysis returned essentially the 
same maximum likelihood (ML) result (unfortunately, surrounded by large remaining uncertainties): the relic gravitational 
waves are present and they are responsible for approximately 20% of the temperature quadrupole. We identify and discuss 
, the reasons by which the contribution of gravitational waves can be overlooked in a data analysis. One of the reasons is a 
. misleading reliance on data from very high multipoles £, another - a too narrow understanding of the problem as the search 
' for B-modes of polarization, rather than the detection of relic gravitational waves with the help of all correlation functions. 
Our analysis of WMAP5 data has led to the identification of a whole family of models characterized by relatively high values 
of the likelihood function. Using the Fisher matrix formalism we formulated forecasts for Planck mission in the context of 
this family of models. We explore in details various 'optimistic', 'pessimistic' and 'dream case' scenarios. We show that in 
some circumstances the B-mode detection may be very inconclusive, at the level of signal-to-noise ratio S/N = 1.75, whereas a 
smarter data analysis can reveal the same gravitational wave signal at S/N = 6.48. The final result is encouraging. Even under 
' q ' ■ unfavourable conditions in terms of instrumental noises and foregrounds, the relic gravitational waves, if they are characterized 
j-^h' by the ML parameters that we found from WMAP5 data, will be detected by Planck at the level S/N = 3.65. 
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I. INTRODUCTION 



A complete cosmological theory is supposed to explain not only the present state of the observed Universe (see, for 
example, Q), but also its early dynamical behaviour and possibly its birth Q. Our present state P is characterized 
by approximate large-scale homogeneity and isotropy within a patch of the size l p « lO 3 ^ « 10 31 cm (see last paper in 
Q) and the averaged energy density p p c 2 of all sorts of matter in this patch p p = 3H 2 /8ttG w 10~ 29 g/cm 3 , where Hq 
is the present day Hubble parameter and Ijj — c/ Hq m 10 28 cm. The limits of applicability of the currently available 
theories are set by the Planck density ppi — c 5 /G 2 h w 10 94 g/cm 3 and the Planck size lp\ — (Gh/c 3 ) 1 / 2 m 10 _33 cm. 
One can imagine that the embryo Universe was created by a quantum-gravity (or by a 'theory-of-everything') process. 
The emerging classical configuration was probably characterized by the near-Planckian energy density and size. The 
total energy, including gravity, was likely to be zero then, and remains zero now. 

The problem is that this hypothesis requires further assumptions. The arising classical configuration can not reach 
the present state P if it expands all the time according to the usual laws of radiation-dominated and matter-dominated 
evolution. By the time the Universe (i.e. the patch of approximate homogeneity and isotropy) has reached the size l p , 
its energy density would have dropped to the level many orders of magnitude lower than the required p p . Therefore, the 
newly born Universe needs a primordial kick before it can join the pathway of normal radiation-dominated expansion. 
The kick should allow the size of the patch to increase by about 33 orders of magnitude without losing too much of 
the energy density of whatever substance that was there, or maybe even slighly increasing this energy density at the 
expense of the energy density of the gravitational field. 

The relic gravitational waves [I[ are necessarily generated by a strong variable gravitational field of the very early 
Universe. They are the cleanest probe of what was happening during the violent times of the initial kick. Specifically, 
the quantum-mechanical Schrodinger evolution transformes the initial no-particle (vacuum) state of the gravitational 
waves into a multiparticle (strongly squeezed vacuum) state. Under certain additional conditions, the same holds 
true for other degrees of freedom of the gravitational field (metric) perturbations, including those representing the 
density perurbations. As a result, the patch of homogeneity and isotropy will necessarily be augmented by primordial 
cosmological perturbations of quantum-mechanical origin. This process is called the superadiabatic, or parametric, 
amplification; for a recent review of the subject, see the last paper in [3]. 

As before (see p} and references there) , we are working with 

ds 2 = -c 2 dt 2 + a 2 (t)(<% + h ij )dx t dx j = a 2 (ri){-dn 2 + (% + h i j)dx i dx j ] 



and the Fourier-expanded gravitational field (metric) perturbations 
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The polarization tensors V%j (n) (s = 1,2) describe either the two transverse-traceless components of gravitational 
waves (gw), or the scalar and longitudinal-longitudinal components of density perturbations (dp). Assuming the 
initial vacuum state of participating perturbations, the resulting metric power spectra h 2 (n,r]) are given by 



IM^ x )^ x )lo) 



— h 2 (n, V ), h 2 (n, V ) = ^n 2 V | h n (r/)| : 

s=l,2 



(2) 



where the mode functions h n (rj) are taken either from gw or dp equations, and C — vT67ripi for gravitational waves 
and C = V247rZpi for density perturbations. 

The numerical levels and shapes of the generated power spectra are determined by the strength and variability 
of the gravitational 'pump' field. The simplest assumption about the initial kick is that its entire duration can be 
described by a single power-law scale factor [J] 



ofa) = lo\v\ l+P , 



(3) 



where l and (3 are constants, [3 < — 1. Then, the generated primordial metric power spectra (for wavelengths longer 
than the Hubble radius at that time) have the universal power-law dependence on the wavenumber n: 



h 2 (n) a n 2 ^+ 2 ). 

It is common to write these metric power spectra separately for gw and dp: 

h 2 (n) (gw) = B 2 n nt , h 2 (n) (dp) = B 2 n n "~ x . 



(4) 
(5) 
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In the case of power-law scale factors (|3|) (or piece- wise power-law scale factors) , the equations for metric perturbations 
representing gravitational waves and density perturbations are exactly the same [5j. Therefore, according to the 
theory of quantum mechanical generation of cosmological perturbations H, the spectral indices are approximately 
equal, n s — 1 = n t = 2((3 + 2), and the amplitudes B t , B s are of the order of magnitude of the ratio Hi/Hp\, where 
Hi ~ c/l a is the characteristic value of the Hubble parameter during the kick. (An initial kick driven by a scalar 
field is usually associated with inflation.) In what follows, we are using the numerical code CAMB Q and related 
notations for gw and dp power spectra adopted there: 

P t (k) = A t (k/k ) n \ P s (k) = A s (k/k Q ) n °-\ (6) 

where k = 0.002MPC" 1 . 

There is no doubt that the metric perturbations with wavelengths greater than the Hubble radius in the times of 
recombination do exist. Indeed, it is known for long time p| that it is precisely this sort of long-wavelength metric 
perturbations that provide the main contribution to the lower multipoles, starting from I = 2, of the Cosmic Microwave 
Background (CMB) temperature anisotropics. The very existence of CMB anisotropies at lower £'s [1, Q testifies to 
the existence of such long- wavelength perturbations. They are likely to be the perturbations of quantum-mechanical 
origin. 

The assumption of a single power-law index (3 is the simplest and easiest to analyze, but it is too strong. It has the 
consequence that one and the same spectral index describes the interval of 30 orders of magnitude of wavelengths in 
the primordial power spectra. In reality, as it appears from our CMB analysis below, even at the span of 2 orders of 
magnitude in terms of wavelengths the spectral index n s is likely to be somewhat different. We will discuss this point 
in more detail in the text of this paper. 

We start (Sec. [D} with a significant improvement of our previous analysis [l[ of the 5-year Wilkinson Microwave 
Anisotropy Probe (WMAP5) TE data at i < 100 flpj|. In contrast to [if, we work directly with both TT and TE 
datasets and impose no restrictions on the (constant) perturbation parameters A t , A s ,n t ,n s , except n s — 1 = n t which 
is implied by the theory of quantum-mechanical generation. We work with the quadrupole ratio R 



C£ T 2 (dp) 



R = r^TT I i \ ( 7 ) 



and the remaining two free parameters A s ,n s . This more general analysis returns essentially the same as before [l[ 
maximum likelihood (ML) values: R — 0.229 and n s — 1.086, with approximately the same as before uncertainties. 
We demonstrate that the data at the multipole I = 2 (dubbed "anomalously low" in the literature) are not to be 
blamed for these determinations. After removal of this data point altogether, the ML values do not change much. 
These improvements and cross-checks make more stable and robust our conclusion 1] that the WMAP5 data do 
contain a hint of presence of relic gravitational waves. 

In Sec. IIIII and Sec. IIVI we show in detail how relic gravitational waves can be overlooked in CMB data analysis. 
In Sec. IIIII we concentrate on one of the reasons, which is the attempt of placing the ever "tighter" constraints on 
gravitational waves by using the data from high ^'s of CMB and large-scale structure surveys. These data have nothing 
to do with gravitational waves and they can mislead the identification of gw contribution. Wc show that even the 
use of CMB data from the adjacent interval 101 < i < 220, where the role of gravitational waves is already small, is 
dangerous. We defer to Sec. IIVI a detailed discussion of another recipe to overlook the relic gravitational waves. This 
is the wide-spread 'obsession' with the detection of -B-modes, rather than the detection of relic gravitational waves 
with the help of all available observational CMB channels. 

The prospects of observing relic gravitational waves by the already deployed Planck satellite [ll[ are analysed in 
great detail in Sec. IIVI We adopt improved evaluations of foregrounds [T3|, [HI], [l4[ and instrumental noises. The 
main thrust of the section is the comparison of the performances of various combinations of observational channels: 
TT + TE + EE + BB, TT + TE + EE, and BB alone. We discuss different models of the foregrounds and their 
subtraction, individual 'optimistic', 'pessimistic' and 'dream case' scenarios, as well as complications in the data 
analysis itself. The final conclusions are formulated not only for the model characterized by the set of ML parameters 
derived from the WMAP5 data, but also for the whole class of models characterized by the high values of the 3- 
parameter likelihood function. We show that there exists plenty of situations where the results from the BB channel 
alone are inconclusive, whereas a smarter data analysis can reveal a significant detection. For other approaches to 
observing relic gravitational waves in the CMB temperature and polarization anisotropies see [IH [l6j . 

The good news is that even under unfavorable conditions, the Planck satellite will see the relic gw signal (assuming 
that it has the WMAP5 maximum likelihood value R = 0.229) at a better than 3<r level. Furthermore, we believe 
that the methods and evaluations of this paper can also be used in ground-based and balloon-borne experiments 
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At the end of the introduction, it is important to stress that the current thinking in this area of science is greatly 
influenced by inflationary understanding of quantum mechanics and relativity: "Quantum fluctuations, usually ob- 
served only on microscopic scales, were stretched to astronomical sizes and promoted to cosmic significance as the 
seeds of large scale structure" [24[ , "the superluminal expansion of space during inflation stretched these scales outside 
of the horizon" (25j . "Inflation... stretched space. ..and promoted microscopic quantum fluctuations to perturbations 
on cosmological scales. Inflation makes detailed predictions..." (26], and so on. 

Indeed, inflationary views on physics have translated into inflationary observational predictions. They are encap- 
sulated in the formula for the predicted scalar metric power spectrum of density perturbations P s ~ Pt/t, which 
is divergent at small e (e = —H/H 2 , e = in the standard inflation), and the detailed 'tensor-to-scalar' ratio r 
(r = Pt/P s ): 

r = 16e = -8n t . (8) 

The widely quoted limits onr,r< 0.22 (95% C.L.) were derived from the likelihood function for r. The analysis 
has resulted in the maximum likelihood value tml = 0. Since P t ss (Hi/Hpi) 2 ^ and tml = 0, one has to decide 
whether the most likely values of density perturbations P s responsible for the data collected and analyzed by the 
WMAP team are infinitely large, or inflationary predictions are wrong. The existing and planned data analyses are 
usually based on the enforced (incorrect) inflationary relation r = — 8n t ; the final physical conclusions are formulated 
in terms of constraints imposed on the (possibly non-existent) scalar field, and so on (see, for example, @, (27|). 

As for the quantity r, there is no doubt that, in general, the inflationary theory can predict for r everything what 
one can possibly ask for (for a review, see [25l| and references there). But the most advanced inflationary theories, 
operating with warped D-brane inflation [2a|. [29]] . D3-brane inflation in warped throats [30], string theory inflation 
[3l| . etc., either "allow a very low tensor amplitude r 10~ 4 ", or lead to the conclusion that "D3-brane inflation 
in Calabi-Yau throats, or in most tori, cannot give rise to an observably-large primordial tensor signal", or to the 
conclusion that r w 10~ 24 and the "existing models of string theory inflation do not predict a detectable level of 
tensor modes" . These conclusions make the search for the inflationary gravitational waves (i.e. relic gravitational 
waves as presented by inflationists) a senseless enterprise. 

Obviously, in this paper, we are not using the inflationary theory and its observational predictions. (These pre- 
dictions are based on the inflationary hat-trick of extracting arbitrarily large scalar metric perturbations P s out of 
vacuum fluctuations of the scalar field. For a more detailed criticism of inflationary theory, see last papers in [ij and 

!•) 



II. IMPROVED EVALUATION OF RELIC GRAVITATIONAL WAVES FROM WMAP TT AND TE DATA 



A. Likelihood functions and summary of the previous results 



Relic gravitational waves compete with density perturbations in generating CMB temperature and polarization 
anisotropics at relatively low multipoles i < 100. For this reason we focus on the WMAP data at £ < 100. As 
before Q, we use the symbols Cj T , Cj E , C EE , Cf B for CMB power spectra and Dj T , Dj E , D EE , Df B for their 
estimators. In this section we ignore the B-mode of polarization because WMAP puts only u ppe r limits on it . 

The variables Dj T , Dj E and Df E obey the Wishart probability density function (pdf) [J, [3^, |H, [I| 



f( D r,Dj E ,D EE ) 



1 \ n/2 n 3 (xy- z 2 )(»- 3 )/ 2 



4(1 - pj){ajaf) 2 J n^T(n/2)T((n I)/2) 

x ex 4" 2 TT 1 7!)((^ + (^"^)}' (9) 

where / s k y is the sky-cut factor, / s k y = 0.85 for WMAP and / s k y = 0.65 for Planck, and n = (21 + I)/ S ky is the 
number of effective degrees of freedom at multipole t. T is the Gamma-function. This pdf contains the variables 
Df x ' (XX' = TT, TE, EE) in quantities x = n(Df T + Nj T ), y = n(D EE + N EE ), z = nDj E . The information on 
the power spectra Cf x is contained in quantities 

CJ E 



-J = \ CJ T + Nf T , of = X C EE + N EE 



Pt 



I(CJ T + NJ T )(C EE + N EE ) 

where Nj T and NF E are the total noise power spectra. 
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We are mostly interested in TT and TE data, so we shall work with the joint pdf for Dj T and Dj E . This pdf is 
derived from ([9]) by integrating over the variable Df E . The resulting pdf has the form 

f(Dj T ,Dj E ) = {2 1 +" 7 rr 2 (|)(l - pl){ajf n {af ) 2 }~* 

1 



X eX P "1 1 3 1 TTTB _ On-tnEXl ~ Of„T\2 I f ' ( 10 ) 



In order to estimate parameters, such as R, A s , n s , from observations one seeks the maximum of the likelihood 
function. The likelihood function is the pdf in which the data (estimates Df x ) are known while the parameters are 
unknown. Up to a normalization constant, the likelihood function is 



C ocH f(Dr,Dj E ) 



for I = 2, • • •, Umax- It can be rewritten as 

" ™ - E{t4 {Jfy + J^f ~ ) +>»((! " A)(^ff)} + & (11) 

where the constant C is chosen to make the maximum value of C equal to 1. 

Our previous analysis [l[ was based on the background ACDM cosmological model as derived in (35|. In addition 
to the relation nt = n s — 1, the perturbation parameters A s ,At were restricted by the observational condition 
1(1 + l)Cj T /2ir\i = io = 840 /iK 2 . One more restriction was supplied by the phenomenological relation n s (R) = 
0.96 + 0.35-R— 0.07i? 2 which indirectly took into account the data on TT anisotropics. The remaining free parameter 
R was subject to the likelihood analysis. This analysis was directly using the 5-year WMAP TE data at multipoles 
2 < £ < 100 The noise power spectra Nf T , N EE were obtained from the information posted at [l(| and were 
presented as graphs in [lj (Fig. 6). 

The likelihood procedure has resulted in R = 0.240±q (68.3% C.L.). For the ML value R = 0.240, the imposed 
restrictions have produced the full set of perturbation parameters: 

R = 0.240, 72 s = 1.040, A s = 2.034 x 10~ 9 , A t = 0.960 x 10~ 9 (12) 

and n t = 0.040. A shift of the parameter R within its confidence interval would automatically produce a change in 
other parameters too. 

In order to avoid any association with inflationary predictions, we are not using the parameter r. However, if r is 
defined as r = A t /A s (definition used by the WMAP team) without implying inflationary formulas ([8]), then one can 
establish a relation between R and r which depends on the background cosmological model and spectral indices f3(| . 
We derived this relation numerically. For a rough comparison of results one can use r « 2R. 



B. Revised analysis of the WMAP5 data 

In this paper, the previous approach [1[ is improved in two main aspects. First, we work directly with both datasets, 
TT and TE. Second, in the likelihood procedure all three parameters R, A s , n s are kept free. (We tried to include nt 
as a free parameter, but this did not change the results except increasing uncertainties around the ML values.) 

As before, the WMAP5 estimates for Dj T ', Dj E at multipoles I < £ max = 100 are taken from [l(|. The noise 
power spectra Nf T , Nf E are the same as derived in [l| (Fig. 6). Numerical evaluations of the CMB power spectra 
are performed with the help of CAMB code [6j . 

The adopted background model is the best-fit ACDM cosmology [9] (ApJS version) with parameters 

n b h 2 = 0.022671^0059. n ch 2 = 0.1131 ± 0.0034, , . 

A = 0.726 ± 0.015, T roion = 0.084 ± 0.016, h = 0.705 ± 0.013. 1 j 

In numerical calculations, we use the central values of these parameters. 

Applying the Markov Chain Monte Carlo (MCMC) method (see, for example, [13,[3^|), we probe the likelihood 
function (| 1 1 1) by 10,000 samples and determine the position of its maximum in 3-dimensional space R, A s ,n s . The 
parameters of our best-fit model, i.e. the maximum likelihood (ML) values of the perturbation parameters, are found 
to be 

R = 0.229, n s = 1.086, A s = 1.920 x 10" 9 (14) 



5 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Likelihood 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Likelihood 



FIG. 1: The projection of 10,000 samples of the 3-dimensional likelihood function onto the R — n 3 (left panel) and R — A s 
(right panel) planes. The black + indicates the parameters listed in (114[l , 

and n t = 0.086. Obviously, these are only the coordinates of the maximum in the parameter space. There are many 
neighbouring points with almost equally large values of the likelihood C. It is difficult to visualize the 3-dimensional 
region around the maximum, so in Fig. [1] we show the projection of the 10,000 sample points on the 2-dimensional 
planes R — n s and R — A s . 

The color of an individual point in Fig. Q] signifies the value of the 3-dimensional likelihood of the corresponding 
sample. The projections of the maximum (|14[) are shown by a black +. The samples with relatively large values of 
the likelihood (red, yellow and green colors) are concentrated along the curve, which projects into relatively straight 
lines (at least, up to R » 0.5): 

n s = 0.98 + 0.46i?, A s = (2.27- 1.53i?) x 10~ 9 . (15) 

These combinations of the parameters R, n s , A s produce roughly equal responses in the CMB power spectra. The 
best-fit model is a particular point on these lines, R — 0.229. We will be using this one-parameter family of 
models (fl"5j) in our study of the detection abilities of the Planck mission in Section IIVI 

Before comparing the new and old results, it is instructive to explore the marginalized 2-dimensional and 1- 
dimensional distributions. The marginalized distribution over a parameter is the integral of the likelihood function 
over that parameter. By integrating C (jlip over A s or n s we derive 2-dimensional likelihoods in R — n s or R — A s 
spaces. Then we apply the standard procedure of finding the maxima and confidence contours. 

In the left panel of Fig. [2] we show the ML point (marked by a red x) and the 68.3% and 95.4% confidence contours 
(red solid lines) in the R — n s plane. The 2-parameter maximum is located at 

R= 0.203, n s = 1.082. (16) 

In the same panel we show 2-dimensional confidence contours as given by the WMAP team @. We transferred their 
contours originally plotted in r — n s plane to our R — n s plane using the numerical relation between R and r. Their 
contours are based either on the assumed constancy of the spectral index n s throughout all the explored multipolcs 
(black dashed curves), or on a simple running of n s , n s (k) — n s (k ) + m (^) (blue dash-dotted curves). 

In the right panel of Fig. [2]we show the ML point (marked by a red x ) and the 68.3% and 95.4% confidence contours 
(red solid lines) in the R — A s plane. The 2-parameter maximum is located at 

R = 0.211, A s = 1.877 x 10~ 9 . (17) 

Integrating the likelihood function L (fTTj) over two parameters (A s , n s ), (A s , R) or (R, n s ), we arrive at 1- 
dimensional distributions for R, n s or A s , respectively. We plot these distributions in Fig. [3] The ML values of these 
parameters and their 68.3% confidence intervals are given by 

R = 0.266 ± 0.171, n s = 1.107±g;g^, A s = (1.768±g;|^) x 10" 9 - (18) 
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FIG. 2: The ML points (red x) and the 68.3% and 95.4% confidence contours (red solid lines) for 2-dimensional likelihoods: 
R — n s (left panel) and R — A 3 (right panel). In the left panel, the WMAP confidence contours are also shown for comparison. 




FIG. 3: 1-dimensional likelihoods for R (left), n 3 (middle) and A s (right). 

Comparing the old, Eq. (fl~2"|) . and new, Eqs. {Tj]) , (fTfjj) . (fTT)) . (fT5|) . results, one can conclude the following. First, 
all the results are close to each other and deviate little around the rigorous 3-dimensional ML values (fT4j) . Second, 
the parameter R persistently indicates a significant amount of relic gravitational waves, even if with a considerable 
uncertainty. The R = hypothesis (no gravitational waves) appears to be away from the R = 0.229 model at about 
la interval, or a little more, but not yet at a significantly larger distance. Third, the spectral indices n s , nt persistently 
point out to the 'blue' shape of the primordial spectra, i.e. n s > l,n t > 0, in the interval of wavelengths responsible 
for the analyzed multipoles I < £ m ax = 100. This puts in doubt the (conventional) scalar fields as a possible driver 
for the initial kick, because the scalar fields cannot support (3 > —2 in Eq. (|3|) and, consequently, n s > l,rtt > in 
Eq. ©. 

C. Quadrupole data and extrapolation of the ML model to higher multipoles 

It is known [391 ] that the actually observed quadrupole Dj]i 2 nas anomalously low value in comparison with the 
usually plotted graphs of the best-fit CMB power spectra. Since our results prefer a somewhat 'blue' primordial 
spectrum, the natural question arises whether the low value of the quadrupole is not the reason entirely responsible for 
our evaluation. In order to answer this question we have conducted the likelihood analysis without using the observed 
data points DjT 2 and Dj^ 2 . We found that even this drastic measure of complete removal of these data points 
does not change our results qualitatively. The parameters of the maximum likelihood model modify to R = 0.190, 
n s = 1.067 and A s = 1.993 x 10~ 9 . This is one more evidence of the stability of indications on the presence of relic 
gravitational waves in the WMAP5 CMB data. 

As was already stressed in the paper, we analyze only those WMAP5 data where one can expect to find relic 
gravitational waves, that is, in the range 2 < t < 100. Therefore, the parameters (fT4|) apply only to wavelengths 
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FIG. 4: The unbinned WMAP5 TT (left panel) and TE (right panel) data [T(J, along with their uncertainties due to noises 
and statistics. The red solid lines show the ML model (|14[) extrapolated with constant spectral indices to high multipoles. For 
comparison, the dashed green line shows the WMAP5 best-fit "no gravitational waves" model. 



responsible for that range. We will show in Sec. IHII that it can be misleading to try to constrain gravitational waves 
by the data outside this interval of multipoles, as the spectral indices may change. Nevertheless, it is interesting to 
see what kind of CMB power spectra the ML model (fT4|) generates, if the spectral indices are assumed fixed at their 
values (|14[) throughout all the relevant wavelengths. 

In Fig. |U we show these extrapolated TT and TE power spectra built on the ML parameters (fT4"]) . One can see 
that these spectra are not too far away from the "no gravitational waves" spectra advocated by the WMAP team. As 
one could expect, the somewhat 'blue' spectral index n s in (|14[) makes the extrapolated spectra positioned somewhat 
above the WMAP5 spectra at very large multipoles. 



III. HOW RELIC GRAVITATIONAL WAVES CAN BE OVERLOOKED IN THE LIKELIHOOD ANALY- 
SIS OF TT AND TE DATA 

With all the reservations already stated, our results are markedly different from the WMAP5 conclusions [9(. The 
WMAP team has found no evidence for gravitational waves and arrived at a red spectral index n s = 0.96. The WMAP 
findings are symbolized by black dashed and blue dash-dotted contours in Fig. O It is important to understand the 
reasons for these disagreements. 

Two differences in data analysis have already been mentioned. We restrict our analysis to multipoles £ < 100, 
whereas the WMAP team uses the data at all multipoles up to £ ~ 1000 keeping spectral indices constant. We use the 
relation n t = n s — 1 implied by the theory of quantum- mechanical generation of cosmological perturbations, whereas 
the WMAP team uses the inflationary 'consistency relation' r = — 8n t which automatically sends r to zero when nt 
approaches zero. There could be some discrepencies in treating the noises, but we think we effectively followed [l| the 
WMAP prescriptions. After several trials, we came to the conclusion (with heavy heart, as Einstein used to say) that 
it is the assumed constancy of spectral indices in a broad spectrum that is mostly responsible for the disagreement, 
and it should be abandoned. The constancy of n s over the vast region of wavenumbers, or possibly a simple running 
of n s , is a usual assumption in a number of works Q, [32l. |40|. 

In order to understand the impact of higher multipole data, we first probed the likelihood function and estimated 
the parameters from the data in the range of multipoles 101 < £ < 220. The procedure was exactly the same as 
was used in the analysis of 2 < £ < 100 data. The maximum of the 3-parameter likelihood function was found at 
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FIG. 5: Same evaluations as in the left panel in Fig. [2] but derived from WMAP5 TT + TE data in the interval 101 < I < 220 
(left panel), and in the interval 2 < I < 220 (right panel). WMAP5 contours [i| are shown for comparison. 




FIG. 6: 1-dimensional likelihood for n s derived from WMAP5 TT + TE data in the interval 101 < £ < 220 (left panel), and in 
the interval 2 < £ < 220 (right panel). 

n s = 0.923, R = 0.022 and A s = 2.65 x 10~ 9 , that is, at a distinctly 'red' spectral index n s . The 2-dimensional 
marginalized distribution R — n s (analogous to the left panel in Fig. [5]) is also different. It is shown in the left panel 
of Fig. El The large uncertainty surrounding R = 0.022 reflects the fact that, at these multipoles, the contribution 
of relic gravitational waves becomes very small. It is the density perturbations that play dominant role here, and at 
higher multipoles. The 1-dimensional marginalized distribution for n s (analogous to the middle panel in Fig. [3]) is 
shown in the left panel of Fig.[6l This distribution gives n s = 0.948^o'q61 (68.3% C.L.). Obviously, this value of n s is 
significantly smaller than the one in Eq. (|18p. and the two evaluations do not overlap in la confidence interval. This 
is a clear indication that the spectral index n s can hardly be treated as one and the same constant throughout all the 
wavelengths responsible for 2 < t < 100 and 101 < I < 220 intervals. 

As one could expect, exactly the same analysis of the whole interval 2 < £ < 220 from the position of constant n s 
leads to the ML values which are intermediate between evaluations at the intervals 2 < I < 100 and 101 < £ < 220 
separately. The 3-parameter likelihood analysis applied to WMAP5 TT + TE data in the interval 2 < £ < 220 has 
resulted in the ML values n s = 0.973, R = 0.019 and A s = 2.39 x 10 -9 . As expected, the ML value n s = 0.973 
is in between the ML values n s — 1.086 and n s = 0.923 from the two adjacent intervals of I. The 2-dimensional 
marginalized constraints are shown in the right panel of Fig. [5j In comparison with the left panel, the uncertainty 
contours are much closer to the WMAP5 evaluations. The 1-dimensional marginalized distribution for n s is shown 
in the right panel of Fig. [5] It gives n s — 0.998^q'q29 (68.3% C.L. ). Again, this is the intermediate value of n s in 
comparison with evaluations from the two intervals of £ separately. Our evaluation of n s from the data in the interval 
2 < £ < 220 still gives a little bit higher value of n s than n s — 0.96 found by the WMAP team, but presumably the 
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remaining difference is accounted for by the multipoles I > 220 and other data sets that were included in the WMAP 
derivation. 

One can see now why the inclusion of data from I > 100 is dangerous. Although these data have nothing to do 
with gravitational waves, they bring n s down. If one assumes that n s is one and the same constant at all £'s, this 
additiosnal 'external' information about n s affects uncertainty about R and brings R down. This is clearly seen, for 
example, in the left panel of Fig. [2] The localization of n s near, say, the line n s = 0.96 would cross the solid red 
contours along that line and would enforce lower, or zero, values of R. However, as we have shown above, the n s is 
sufficiently different even at the span of the two neighbouring intervals of ^'s that we discussed. These considerations 
as for how relic gravitational waves can be overlooked have general significance and will be applicable to any CMB 
data and data analysis. 

IV. FORECASTS FOR THE PLANCK MISSION 

Now, that Planck satellite [ill ] is launched, it is important to see in detail how the maximum likelihood model 
(P3J found from WMAP5 data, as well as neighbouring models defined by Eq. (fTS)) . will fare in the less noisy Planck 
data. In our previous work [l[, we restricted our consideration to only two separate information channels, namely to 
TE and BB correlation functions. We took into account only one frequency channel 143GHz and its instrumental 
noise. We replaced the unknown level of residual foregrounds and systematics by the increased total noise in the BB 
channel, calling it a 'realistic' BB channel. On the data analysis side, we performed the likelihood analysis with a 
single parameter R, assuming that other parameters A s , n s and n t = n s — 1 are known, as long as they are related 
to R by the imposed restrictions. 

In this paper, we study the detection abilities of Planck in a much more comprehensive manner. First, we consider 
all available information channels, i.e. TT, TE, EE and BB correlation functions, and their combinations. Second, 
in evaluating the instrumental noises, we take into account all the Planck's frequency channels (either three, or even 
all seven). Third, we include in the total noise the residual foreground contamination from the synchrotron and 
dust emissions, and introduce the 'pessimistic' case when the foregrounds are not removed at all while the nominal 
instrumental noise in the BB channel is increased by a factor of 4. Finally, we search not only for a single parameter 
R, but where the computational resources allow, we evaluate the uncertainties for R arising in the procedure of joint 
determination of all three unknown parameters R, n s , A s from a given set of observational data. 

A. The signal-to-noise ratio in the measurement of R 

1. The noise power spectra and the definition of S/N 

Surely, in the focus of our attention is the detection of relic gravitational waves. Their CMB contribution is 
parameterized by R. To quantify the detection ability of a CMB experiment, we introduce the signal-to-noise ratio 
1 

S/N = R/AR, (19) 

where the numerator is the true value of the parameter R (or its ML value, or the input value in a numerical simulation) 
while AR in the denominator is the expected uncertainty in determination of R from the data. 

In formulating the observational forecasts, it is common to use the Fisher matrix formalism. We outline its main 
features in Appendix [Al The uncertainty AR, for a given R, depends on noises, statistics of the searched for CMB 
signal (which is random in itself), and the number of parameters subject to evaluation from a given dataset. The 
instrumental and foreground noises comprise the total effective noise power spectra N xx which enter the Fisher 
matrix (|A5[) and its element Fan, Eq. (IA7|) . as well as the element RR of the inverse matrix. Depending on the 

number of sought after parameters, one calculates AR either according to Eq. (|A6|) or according to Eq. (|A8|) . 

The noise power spectra N xx are explained in Appendix [Bj We ignore the cross-correlated noises, i.e. when 
X 7^ X' , so we are working only with N xx . The instrumental characteristics of all seven frequency channels [ll[ (we 
mark them by symbol i = 1, • • •, 7) are listed in Table [J of Appendix [Bl where LFI means Low Frequency Instrument 
and HFI High Frequency Instrument. Following |ll[ we treat only three frequency channels 100GHz, 143GHz and 
217GHz as providing CMB data, whereas other channels are supposed to be used for determining the foregrounds. 
However, in Section IIV CI we study the improvements that would arise if one could use all seven frequency channels 
for CMB analysis. It is seen from (|Bip . that three channels are better than one, and seven channels are better than 
three. 
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FIG. 7: The total noise power spectrum £(£ + l)N BB /2iv in comparison with the BB power spectrum 1(1 + l)Cf B /2ir for the 
model 03) with R = 0.229. 

The residual foreground noise N xx {i) adds to the instrumental noise N^^(i) and increases the total noise N xx 
(|B1I) . We neglect the effect of foregrounds in TT channel because this noise is small in comparison with the signal [4JJ . 
The foreground contamination cannot be neglected in polarization channels. The synchrotron and dust emissions are 
expected to be the dominant hindrances in the Planck frequency range [HI, HH, |43| . The charecteristics of the accepted 
foreground models are listed in Appendix [BJ In what follows, we are using the more severe Dust A model, whereas 
the more favorable Dust B model is considered only in Section llV CI 

The ways of mitigating the foreground contamination are discussed in a number of papers 

ma mm gam. 

In this paper, we take a phenomenological approach and quantify the residual noise by the factor er fg (see }13| and 
Appendix|B]) which multiplies the model power spectra C x f (i), C xx (i) of the synchrotron (S) and dust (D) emissions. 
We consider three cases, <r fg = 1 (no foreground removal), er fg = 0.1 (10% foreground residual noise) and cr fg = 0.01 
(1% foreground residual noise). In order to gauge the worst case scenario, we consider also the 'pessimistic' case, 
where <r fg = 1 and NP B t (%) at each frequency v% is 4 times larger than the values listed in Table Q] This increased 
noise is meant to mimic the situation where it is not possible to get rid of various systematic effects (49|, the E-B 
mixture cosmic lensing [5l|, etc. which all affect the BB channel. 

To illustrate the expected total noise, including the different levels of foreground contamination, we show in Fig. [7] 
the total noise power spectrum N BB calculated according to Eq. (|Bl[) . For reference, the black curve shows the power 
spectrum Cf B for the maximum likelihood model with R — 0.229 (see Eq. (fl"4"])) extrapolated up to £ = 1000. It 
is seen from the graphs that the role of foregrounds is restricted to relatively small multipolcs £ < 20. For higher 
multipoles, the total noise N BB is dominated by the instrumental noise and does not depend on cr fg . It can also 
be seen from Fig. [7] that it is only for small values of <7 fg , i.e. in the case when the foreground contamination is 
significantly suppressed, that the signal Cf B will be greater than the noise N BB at lowest multipoles. Thus, even 
in the case of small <7 fg , the Planck's BB channel will be mostly sensitive to gravitational waves in the epoch of 
rcionization. For larger values of cr fg , the relative contribution of lower multipoles to the total S/N is diminished, 
signifying the overall reduction of S/N . In this case, the main contribution to S/N comes from higher multipolcs 
I > 20 and S/N will ultimately be limited by the level of instrumental noise. 
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FIG. 8: The S/N as a function of R. In the left panel, the uncertainty A_R is calculated according to Eq. (|A6|I . whereas in the 
right panel it is calculated according to Eq. (|A8[) . 



2. The analysis of the family of models S15}) . including the ML model \1J$ 

The set of maximum likelihood parameters (|14| is the best set among the 'almost equally good' sets, defined by 
([15]) . In a sense, Eq. (fT5|) is our theoretical prediction, based on the analysis of WMAP5 data, of the best viable 
perturbation models. This family of models is the subject of the Fisher matrix analysis below. With all the noises 
N xx defined by (|Blj) and all the power spectra C xx calculable from the family parameters (fT5|) . we have enough 
ammunition to find the quantities (IA5j) . 

We start from the idealized situation, where only one parameter R is considered unknown and therefore the 
uncertainty AR can be found from (|A6[) . All the information channels, TT, TE, EE, BB, are used in the calculation 
of Frr, Eq. (|A7|) . The results for S/N are plotted in the left panel of Fig.[H Four options are depicted, a = 0.01, 0.1, 1 
and the pessimistic case. 

The results for the benchmark model (TTi|) are given by the intersection points along the vertical line R — 0.229. 
The signal to noise is high, S/N — 9.35, if the foregrounds can be removed to the level cr fg = 0.01, and S/N decreases 
to S/N = 8.80, S/N = 7.72, S/N = 6.48 for <j fg = 0.1, er fg = 1, and the pessimistic case, respectively. In all these 
cases, the S/N is impressively large. A value of R down to R — 0.03 can be measured with S/N = 2 in the optimistic 
scenario er fg = 0.01, whereas S/N > 2 is achieved for R > 0.064 in the pessimistic case. Typically, the optimistic 
scenario gives S/N about 1.5 times larger than the pessimistic case, with the disparity growing larger for smaller 
values of R. 

As one could expect, the uncertainty AR grows and S/N decreases in the realistic situation, when all unknown 
parameters R, n s and A s are supposed to be evaluated from one and the same dataset. In this case, AR should be 
calculated according to (IA8|) . Again, calculating RR , we take into account all the information channels TT, 

TE, EE, BB. The results for S/N are presented in the right panel of Fig. [H For the benchmark model (fl4|) . the 
signal-to-noise ratios are smaller than in the left panel: S/N — 6.69, 6.20, 5.15 if <7 fg = 0.01, 0.1, 1, respectively. 

However, the good news is that even in the pessimistic case one gets S/N > 2 for R > 0.11, and the Planck satellite 
will be capable of seeing the ML signal R = 0.229 at the level S/N — 3.65. 
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B. Contributions of individual information channels and individual multipoles to the total S/N 



The evaluations in Sec. IIV A~2l assume that all the correlation functions, TT+TE + EE+BB, are taken into account 
and all the relevant multipoles £ participate in the summations. The total S/N can only be worse if something is 
missing either in the information channels or in the accessible multipoles. In order to gain further insight into the 
detection ability of Planck, it is instructive to make a break-up of S/N over combinations of information channels 
and multipoles. It is unclear how to do this in the case of Eq. (|A8|) . but it is relatively easy to do this in the case of 
Eq. (|A6|) . We will limit ourselves to this latter case which is sufficient for the purpose of illustration. 

The S/N from Eq. (fT!))) , together with AR from Eq. (|A6[) . can be rewritten as 

(S/N) 2 = R 2 F RR , (20) 
The F RR from Eq. (|A7[) is a simple sum over multipoles £, so the (S/N) 2 in ((20| can be decomposed into individul 



-contributions 



(S/N) 2 = J2(S/N) 2 . (21) 



t 



Each individul (S/N) 2 depends on all information channels, with the _Bi?-channel factored out, as was infixed from 
the very beginning in the form of Eq. (|A2j) , 



(Q/AT\2 — R 2 (2£ + l)/sky f tn TT\2 <^Cf E dCf E <r EE\1^t 9Cj T 



2 (Cj T Cf E - (Cj E ) 2 ) I dR dR OR dR 



dCj E dCj E . ,„ TE ,2 dCj T dC EE 



(pTTnEE 1 lnTE\ 2 \ ^Lg , 



+2 

+4Ll {Ce dR +Ce dR ) dR J 2(Cf B ) 2 9i? 9i? ' 



dR dR dR dR 

,TE( r T T dCf E , r EE dC J T \ dCf E \ , i? 2 (2£+l)/ sky aCf s 9Cr 



The natural break-up of Eq. (J22J) into combinations of the information channels is TT+TE+EE+BB, TT+TE+EE 
and _B£? alone. 



Signal to noise ratio for different combinations of information channels 

Using either all terms in Eq. ([22|) . or everything without BB, or BB alone, we calculate (S/N) 2 , Eq. (f2Tj) . for three 
combinations of channels: TT + TE + EE + BB, TT + TE + EE and BB alone. The (S/N) 2 for the first (full) 
combination is the sum of (S/N) 2 for the other two. Since the foreground removal is a major concern, we separate 
the results into four groups - <7 fg = 0.01, 0.1, 1 and the pessimistic case. The results for S/N are shown in Fig. [H] in 
four panels. Certainly, the (red) lines marked by TT + TE + EE + BB in four panels are the same lines that are 
shown in the left panel of Fig. [5] for the corresponding case. The copy of the line representing the optimistic scenario 
in the upper left panel, i.e. TT + TE + EE + BB together with cr fg = 0.01, is shown by a (black) dashed line in other 
panels as a reminder of what can be achieved in the optimistic scenario. 

Surely, the combination TT + TE + EE + BB is more sensitive than any of the other two, TT + TE + EE and 
BB alone. For example, in the case cr fg = 0.1, the use of all channels provides S/N which is ~ 50% greater than BB 
alone and ~ 30% greater than TT + TE + EE. The situation is even more peculiar in the pessimistic case. The ML 
model (|14p can be barely seen through the popular B-modes alone, because the BB channel gives only S/N = 1.75, 
whereas the use of all channels can provide a confident detection with S/N = 6.48. 

Comparing TT + TE + EE with BB, one can see that the first method is better than the second, except in the 
case when <7 fg = 0.01 and R is small (R < 0.16). In the pessimistic case, the role of the BB channel is so small that 
the TT + TE + EE method provides essentially the same sensitivity as all channels TT + TE + EE + BB together. 

Considering the S/N for BB alone, it is worth noting that since the BB channel is not sensitive to A s and n s 
the uncertainty (|A8[) reduces to (|A6[) . Therefore, although the results for BB channel alone, shown in Fig. [SI were 
derived under the assumption of a single unknown parameter R, they apply also in the general case when all three 
parameters, R, n s , A s , are supposed to be determined from the temperature and polarization data. 



2. Multi-pole decomposition of S/N 

It is seen from Eq. (f2"Tj) the the total (S/N) 2 is a sum over ^-contributions (S/N) 2 given by Eq. (f2"2")l . Formally, the 
sum can extend to large fs, but a smaller i?-signal and larger noises make the large ts irrelevant anyway. We go up 
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FIG. 9: The S/N for different combinations of the information channels, TT + TE + EE + BB, TT + TE + EE and BB. 



tO (.max = 100. 

In Fig. [TO] and Fig. [11] we plot the quantity (S/N) 2 as a function of £ for different combinations of information 
channels and different er fg , including the pessimistic case. Fig. [TO] presents the multipole decomposition for the ML 
model ([14| with R = 0.229, while Fig. E] sorts out the R = 0.05 model characterized by S/N = 3 for cr fg = 0.01 (see 
left panel in Fig. [8]). Note the differences in scaling on vertical axes in Fig. [TO] and Fig. [TT] 

Both figures show again that a very deep foreground cleaning, a 6 = 0.01, makes the very low (reionization) 
multipoles ( ~ 10 the major contributors to the total (S/N) 2 , and mostly from the BB channel. This is especially 
true for the lower-i? model R = 0.05. However, for large er fg = 0.1, 1, and especially in the pessimistic case (see 
the lower right panels in Figs. [TO1 and ITT]) . the role of the BB channel becomes very small at all £'s. These detailed 
illustrations in terms of multipole decomposition of (S/N) 2 are of course fully consistent with the integrated results 
of Sec. IIVB1I 

At the same time, as Figs. [TO] and [TT] illustrate, the ^-decomposition of (S/N) 2 for TT + TE + EE combination 
depends only weakly on the level of cr fg . Furthermore, in this method, the signal-to-noise curves generally peak at 
I ~ (20 — 50). Therefore, it will be particularly important for Planck mission to avoid excessive noises in this region 
of multipoles. 
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FIG. 10: The individual (S/N)J as functions of £ for different combinations of information channels and different levels of 
foreground contamination. Calculations are done for the ML model (|14[) with R = 0.229. 



C. On the possibilities to get a better value for S/N 

So far, the S/N was calculated using the generally quoted realistic assumptions about instrumental and environmen- 
tal noises. It appears that some reserves to get better values for S/N still exist. These possibilities seem speculative 
but worth exploring. First, one may find a way of using the outputs of all seven frequency channels (listed in Table HJ 
for the CMB analysis. The summation over i = 1, • • •, 7, instead of i = 1, 2, 3, would effectively reduce the instrumental 
part of noise in N xx (|B1I) . Second, we may be lucky and the Dust B model turns out to be correct one, rather than 
the more severe Dust A model. This would reduce the foreground part of noise in N xx (|B1|) . 

Below, we consider three possibilities of these improvements in S/N . Namely, the use of 3 frequency channels and 
the validity of Dust B model, the use of 7 frequency channels and the validity of Dust A model, and the 'dream case' 
of using all 7 frequency channels in the conditions of validity of Dust B model. For these three options, we do exactly 
the same calculations that were done in Sec. II V A~2l and depicted in Fig. [8l The left panel in Fig. [8] translates into the 
three upper panels in Fig. [T^l and the right panel in Fig. [3] translates into the three lower panels in Fig. [T2l 

Certainly, one sees considerable improvements in S/N, especially in the 'dream case' scenario (upper right and 
lower right panels in Fig. I12[) . For example, concentrating on the solid line in the lower right panel, one finds that the 
ML model dHJ) with R = 0.229 would be detectable at the level S/N = 5.11, instead of S/N = 3.65 that we found in 
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FIG. 11: The same illustrations as in Fig. OH but for the model JEJ with R = 0.05. 

the right panel in Fig. [8l This would significantly increase the chance of observing relic gravitational waves with the 
help of the Planck mission. 

V. CONCLUSIONS 

Being in the possession of general theoretical confidence that the relic gravitational waves are expected to be present 
in the observed CMB anisotropics, we improved our previous analysis of WMAP5 data and made comprehensive 
forecasts for the Planck mission. 

The improvements in the analysis of the WMAP5 TT and TE observations made our approach more general and 
stable. The new analysis avoids phenomenological relations and restrictions, and deals directly with the complete 
3-parameter likelihood function and its marginalizations. The result of this analysis is very close to the previous 
one: the maximum likelihood value for the quadrupole ratio R is R = 0.229. This means that approximately 20% 
of the temperature quadrupole is caused by gravitational waves and 80% by density perturbations. Although the 
uncertainties due to large noises are still large, this result can be regarded as an indication of the presence of relic 
gravitational waves in the lower I CMB anisotropies (we would love to call it a suspected detection, but we resist this 
temptation). 
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FIG. 12: The improved values of S/N. The three upper panels should be compared with the left panel in Fig. [8] and the three 
lower panels with the right panel in Fig. [8] 



We identify and study in detail the reasons by which the contribution of relic gravitational waves can be overlooked 
in a data analysis. One of the reasons is the unjustified reliance on constancy of the spectral index n s and the inclusion 
of data from very high multipoles £. Another reason - a too narrow understanding of the problem as the search for 
B-modes of polarization, rather than the detection of relic gravitational wave with the help of all correlation functions. 

Our forecasts for Planck are also based on our analysis of WMAP5 observations. We identify the whole family 
of models, that is, the sets of perturbation parameters R, n s , A s , which are almost as good as the unique model 
with the maximum likelihood, C = 1, set of parameters R = 0.229, n s = 1.086, A s = 1.920 x 10~ 9 . For the same 
observational data, these sets of parameters return reasonably high numerical values of the likelihood function C. We 
formulate our forecasts for this family of models, which is characterized by observationally preferred sets of parameters, 
rather than choosing the models and parameters at random and blindly. Our forecasts, based on the Fisher matrix 
techniques, refer to achievable signal-to-noise ratios S/N. We analyze many sources of noise and explore various 
'optimistic', 'pessimistic' and 'dream case' scenarios. We discuss the circumstances in which the B-mode detection is 
very inconclusive, at the level S/N = 1.75, whereas a smarter data analysis reveals the same gravitational wave signal 
at S/N = 6.48. 

The very encouraging final result is that, even under unfavourable conditions in terms of instrumental noises 
and foregrounds, the relic gravitational waves, if they are characterized by the WMAP5 maximum likelihood value 
R = 0.229, will be detected by Planck at the level S/N = 3.65. 
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APPENDIX A: FISHER INFORMATION MATRIX IN THE CMB ANALYSIS 



In the CMB data analysis, the general form of the likelihood function, up to a normalization constant, is 



(Al) 



Since Df is independent of the rest of variables Dj , Df , Df the likelihood function factorizes, 



Cot\[f{Dj T ,Dj E ,D EE )f{D BB ). 



(A2) 



The pdf f(Dj T , Dj E ,Df E ) is the Wishart distribution given by Eq. 



lV ( n -2)/2 e -V/2 



2 n / 2 T(n/2)(C BB 



+ N BB y 



whereas f(Df B ) is the \ 2 distribution: 

(A3) 



where V is V = n(Df B + N BB ) / \Cf B + N BB ). The likelihood function is a function of the sought after parameters 
Si, which in our case are perturbation parameters A t , A s , n t , n s . They enter the likelihood function through their 
presence in the power spectra C xx . The BB part of C depends only on A t and n t . We have reduced the set of 
parameters s,; to R, n s , A s . 

The Fisher information matrix is a measure of the width and shape of the likelihood function, as a function of the 
parameters, around its maximum. The Fisher matrix formalism is used for estimation of the accuracy with which the 
parameters of interest can be found from the data [12, HH . The elements of the matrix are expectation values of the 
second derivatives of logarithm of the likelihood function with respect to s,; , 



d 2 (ln£) 

dsids-j 



(A4) 



where Si are the true values of S; (i.e. values where the average of the first derivative of ln£ vanishes). The angle 
brackets denote the integration over the joint pdf for all Df x . 
Inserting (|A2[) into (|A4[) , one can calculate the Fisher matrix, 
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(A5) 



where Cov 1 is the inverse of the covariance matrix. (Result (|A5|) coincides with Eq. (7) in [53{.) The non- vanishing 
components of the covariance matrix are given by 
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When a particular parameter s is estimated from the data, while other parameters are assumed to be known, the 
uncertainty in the determination of this parameter is given by As = 1/V F ss . However, if all parameters are estimated 
from the data, the uncertainty rises to As = yj (i 7 ' _1 ) ss . The second uncertainty is always larger than the first one or 
equal to it. 

In this work, we are mostly interested in the parameter R, which quantifies the contribution of relic gravitational 
waves to the CMB anisotropics. The definition of the signal-to-noise ratio S/N in Eq. (fT9j) requires the specification 
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of the uncertainty AR. If other parameters are fixed and only R is derived from the data, this uncertainty is given 
by the matrix element Frr, 



AR = 1/VF RR . (A6) 
Explicit expression for Frr follows from Eq. (|A5[) . 
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where Cf x = C xx + N xx , and Cf E = Cj E . 

If other parameters are being determined together with R, the uncertainty AR is given by the RR element of the 
inverse matrix, 



AR 



yJiF-^RR- (A8) 



The uncertainty (|A8|) is always greater than (|A6|) or equal to it. We use (|A6|) and (|A8|) for numerical evaluation of 
S/N. 



APPENDIX B: THE INSTRUMENTAL AND ENVIRONMENTAL NOISE POWER SPECTRA 

CMB experiments use several frequency channels (distinguished by the label i) which have specific instrumental 
and environmental noises at each frequency v^. The optimal combination of the channels gives the total effective noise 
power spectrum Nf x (I3L [iH IH3| 
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where N^ x e (i) and N^ x (i) are the instrumental and residual foreground noise power spectra, respectively. The total 
effective noise power spectrum N xx enters the Fisher matrix (|A5|) . 

In the evaluation of noise power spectra we set the window function equal to 1, which is a good approximation for the 
multipoles considered in this paper, t < 220. The instrumental characteristics of the Planck's frequency channels are 
reported in Table|T]based on [llj. This Table provides N^ x e (i). As for the polarized foregrounds (XX = EE, BB), we 
focus on the synchrotron (S) and dust (D) emissions. The foreground contamination is quantified by the parameter 
<7 fg = 1, 0.1, 0.01 which multiplies the power spectra C x f(i), C E x (i) of the accepted foreground models. The 
smaller cr fg the deeper cleaning. The residual foreground noise is given by [l3j 



*£f(0= E [cf?(i)^+Mff(i)) 



(B2) 



f=S,D 



where Af xx (i) is the noise power spectrum arising from the cleaning procedure itself in the presence of instrumental 



noise. 



Following [13, UjL |42| , for the I and Vi dependences of the synchrotron and dust emissions we take 



^0 



/3s 



(B3) 



and 



2a 1 



e hv /kT _ 2 



ohm/KT _ I 



(B4) 
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TABLE I: Summary of Planck Instrumental Characteristics 



Instrument Characteristic 


LFI 


HFI 


Center Frequency [GHz] 


30 


44 


70 


100 


143 


217 


353 


Angular Resolution [FWHM arcmin] 


33 


24 


14 


10 


7.1 


5.0 


5.0 


[lCr 


27.37 


26.38 


27.20 


3.93 


1.53 


3.62 


33.94 


Ngft and Ngft [10"VK^] 


53.65 


55.05 


55.28 


10.05 


5.58 


15.09 


139.50 



TABLE II: Assumptions about foregrounds [l4| 



Parameter 


Synchrotron 


Dust A 


Dust B 


A S ,D 


4.7 x 10~ b [iK J 


1.0 pK 2 - 


1.2 x 10~ 4 nK J 


v 


30 GHz 


94 GHz 


94 GHz 


to 


350 


10 


900 


a 


-3 


2.2 


2.2 




-2.6 


-2.5 


-1.3 


(3 HH 


-2.6 


-2.5 


-1.4 



In Eq. (|B4|) . p is the dust polarization fraction, estimated to be 5% [13| . and T is the temperature of the dust grains, 
assumed to be constant across the sky T = 18K [Hj]. Other parameters in Eqs. ([Eg]) . (|B4|) are specified in Table HQ 
taken from 

The noise term Nff(i) (/ = S,D) entering Eq. jB2| was calculated in [HE3 



XXf-\ / \ 2a 



n c han(«chan - lj/4 V^ref 



Here, n c h a n is the total number of frequency channels used in making the foreground template map, and v TC f is the 
frequency of the reference channel. In the case of the dust, is rc f is the highest frequency channel included in the 
template making, while in the case of the synchrotron, v rc { is the lowest frequency channel. The value of a is given 
in Table QT] for different foreground models. 

All components of noise are used in numerical calculations of the total noise (|B1[) . 
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